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Abstract 

It is shown that in the physically interesting domain of the parameter space of SU(5) 
supergravity GUT, the Higgs and the Z poles dominate the LSP annihilation. Here 
the naive analyses on thermal averaging breaks down and formulae are derived which 
give a rigorous treatment over the poles. These results are then used to show that 
there exist significant domains in the parameter space where the constraints of proton 
stability and cosmology are simultaneously satisfied. New upper limits on light particle 
masses are obtained. 
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I. INTRODUCTION : Recently there have been extensive investigations of SU(5) 
supergravity models^"^ with electro-weak symmetry broken via radiative corrections^'^. 

Analyses of Refs (6)-(7) were carried out in the framework of No-Scale models^ while the 
analysis of Rcfs (1-5) are for the standard SU(5) supergravity case^*^. In this letter we 
shall discuss only the standard SU(5) Model^^'^. Here, after fixing the Z-boson mass the 
model depends on four arbitrary parameters, aside from the top quark mass rrit, which 
may be chosen to be rrio (the universal scalar mass), mi/2 (the universal gaugino mass), 
Ao (the cubic soft SUSY breaking parameter) and tan P =< H2 > / < Hi > where 
< Hi > gives mass to the down quarks and the leptons and < H2 > gives mass to the 
up quarks. The analyses of Refs (1-4) investigated the full parameter space of the theory, 
Ref (5) investigated the space under one more constraint {Bq = ^0 — ™o where Bq is the 
quadratic soft SUSY breaking parameter) while Refs (1-3) also included in the analysis 
the constraint of proton stability^^. The inclusion of proton stability constraints were 
seen to lead to a number of simple mass relations among the neutralino, chargino and 
gluino mass spectra^"^. One finds for most of the parameter space 2mz^ = m^j = m^^, 
and ~ (l/4)mg for > and ~ (l/3)mg for < 0. (Here Zi^2 are the two 

lightest neutralinos, Wi is the lightest chargino and g is the gluino.) Thus the gluino mass 
(approximately) determines the light neutralino and chargino spectrum. 

Remarkably the standard SU(5) supergravity model under the constraint of proton 
stability also leads to the prediction that the lightest neutralino is the lightest supersym- 
metric particle (LSP)^~^. We investigate here the implications on the parameter space of 
the constraint that the relic density of the lightest neutralino not overdose the universe 
i.e. n^^/i^ < 1, where O^^ = p^^/pc, with the matter density of the lightest neutralino 
Zi and pc = 3HQ/{87rGN) the critical density. Here Hq = hx 100km/ (s.Mpc) and h is the 
Hubble parameter with | < h < 1. Recently it was pointed out that the dominant anni- 
hilation of neutralinos occurs near the lightest neutral Higgs (h) pole in the s-channel for 
the domain of the parameter space satisfying CDF, LEP and proton stability constraints 
and the finetuning requirement that mq^g < ITeV. It is known that the expansion in powers 
of V of the thermally averaged quantity < av > (where a is the spin averaged annihilation 
cross-section of two neutralinos, v is the relative velocity defined by ^/s = 2mz^ + ^ms^v'^, 
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and s is square of the center-of-mass energy), breaks down when ^/s is in the vicinity of 
a pole-*^^. In this case, a careful treatment of integration over the pole in the annihilation 
channel is needed. However, the rigorous analysis even in the non-relativistic approxima- 
tion involves a double integration over the pole (which is numerically intricate) for the 
quantity J = J^^ dx < av > needed to calculate the relic density, {xf = kTf/rrizj^, where 
Tf is the freeze out temperature.) Here we derive rigorous formulas where the integra- 
tions over one of the variables is analytically carried out and the remaining integration is 
smooth over the pole. The analysis here is complete and includes the direct channel Higgs 
and Z-poles as well as t-channel fermion exchange diagrams. Using the rigorous analysis 
for J outlined above we explore the full five dimensional parameter space of the theory 
characterized by mo, mi/2, tan P and rrit under the combined constraints of CDF and 
LEP data, proton stability and relic density. We show that while the parameter space 
is strongly constrained, significant domains in the parameter space remain where all the 
constraints mentioned above are satisfied. Also new limits on the light Higgs, the light 
chargino and on the LSP result. 

II. BASIC FORMULAE : We follow standard procedure^'^ and write the equation 
governing the number density n at time t of the lightest neutralino Zi in a Friedman- 
Robertson- Walker universe with isotropic mass density in the form 

where / = n/T^, x = kT/m^^ {k is the Boltzman constant), Np is the number of degrees 
of freedom at temperature T, G^v is the Newtonian constant and /o = tlq/T^ where uq is 
the number density at thermal equilibrium. The relic density of the LSP is then given by 
the following (approximate) formula^^: 

p,. = 4.75 X 10-« f ' ' r,, (^) J> (2) 

'^^^ J \2.75°K J \J{xf)Jcm^ 

where (T^^/Xy)^ is a reheating factor, is the current temperature and J{xf) is given by 
J{xf) = Jq-'^ < av > dx and: 

(>00 r-OO 

<av>= J dvv\av)e-'''/^y J dw^e-"^'^'' (3) 
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The freezeout temperature Tf is determined by the relation 
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In 



x: 



} <av> V45m5i /(47r'*iV^G'^) 



(4) 



In Eq. (4) < av > is the thermally average of av evaluated at Xf. 

III. INTEGRATION OVER HIGGS AND Z-POLES: J{xf) appearing in Eq.(2) 
can be decomposed as J = Jh+Jz+Jsf where Jh, Jz are the contributions of the s-channel 
Higgs and Z poles, and Jsf is the t-channel contribution from the exchange of squarks and 
sleptons. In the domain of physical interest with finetuning constraints Trtq^g < ITeV, only 
the lightest neutral Higgs h makes a significant contribution to the cross-section. For the 
Higgs pole, using the non-relativistic approximation, we write av in the form 



In Eq.(5) eh = (m^ — 4mU/m|j and 7/^ = mhTh/m1^ where nih is the Higgs mass and 
Th is the Higgs decay width and Ah is^^ 
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Here sin2Q; = — (m^ + m|)(m^ — mf^)^^ sin 2/3 and where uia is the mass of the CP-odd 
Higgs and nin is the mass of the CP-even heavy neutral Higgs. Ci is a color factor which 
is (3,1) for (quarks, leptons) and nu are components of Zi eigen-vector in the basis defined 
in Ellis et al in Ref (14). Using Eq.(5) one can carry out the ^-integration in Jh{xf) and 
get 



Jh{Xf) = 



At 



2V2m^ 



21 



Ilh + —hh 



(7) 



where 
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A similar analysis can be carried out for the Z-Pole and here one finds 
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where Iiz and I2Z are defined anologously to Iih and l2h with m/j, replaced by M^, V z- 
In Eq.(lO) Az IS given by 
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where we have retained only the dominant 6, c and r-contributions, while Bz (in the 
zero-fermion mass approximation) is 



Bz = 



^ ^2 2(2 2 \2 
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In the vicinity of the Higgs (or Z-Pole), J^/ is typically much smaller than J/j (or Jz^ and 
thus we shall use the conventional approximation^^" of axj + \bx^^ in computing Jg/. 
IV. ANALYSIS AND RESULTS : We begin by exhibiting the result that the com- 
putation of Japprox using powcr expansion in v'^ on < crv > is a poor approximation to the 
full analysis of J where rigorous thermal averaging on the Higgs and Z-Poles is carried out. 



The ratio of ^approx/^ — J I Ji 



' approx 



is exhibited in Fig. 1. The results of Fig. 1 show that 



approx 



can be inaccurate by up to 3 orders of magnitude and show a total breakdown of 



the approximate result near the Higgs pole or Z pole. 

To proceed further we must include proton stability constraints. In supergravity 
SU(5), the dominant proton decay proceeds via dimension five operators and involves the 
Higgs color triplet exchange. The most dominant decay mode is p ^ VK^ , and proton 
stability may be conveniently characterized by the value of the dressing loop function B 
that enters in p — > uK'^ decay and is defined in Ref 1. The current Kamiokande bound 



of^^ t{p vK^) > 1 X lO^^yr translates to a bound on B of^''' 

^ < (^) ^"^^ 

where Mh^ is the Higgs triplet mass and Mq is the GUT mass. We also note that the 
simplest GUT sector in SU(5)-'^® leads to the relation Mn^/My = (ax/ao)^ between the 
Higgs triplet mass M^^ and the massive vector boson mass My- [Here a\ = A|/47r and 
A2 enters the GUT superpotential via the term A2-ffi(S + 3M)H2 where Hi, H2 are the 
5, 5 and S is the 24-plet representation of SU(5).] An upper limit on the Higgs triplet 
mass emerges if one assumes that the Yukawa couplings be perturbative at the GUT 
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scale. Estimates on M^g that lead to perturbative A2 lie in the range M^g < 3Mq' to 
Mffg < lOM^. Here as a guideline we shall use a benchmark limit of M^g < GMq- 

We discuss now the result of the analysis. We start at the GUT scale with SU(5) 
supergravity boundary conditions and use the renormalization group equations to evolve 
masses and coupling constants to low energy where a radiative breaking of the electro- 
weak symmetry is achieved. Solutions are subjected to constraints of the CDF and LEP 
data which give lower limits on the SUSY mass spectra, the proton decay constraint of 
Eq.(13) with Mh.^ < GMq, and the relic density constraints discussed in sees I-III. We 
shall also impose the fine tuning condition ruq^g < ITeV. The analysis shows that there 
exists significant domains in the parameter space for both /i > and /x < (/i is the 
Higgs mixing parameter which enters the superpotential via the term ^HiH2) where all 
the desired constraints are satisfied. The allowed parameter space is found to be larger for 
the case n > 0. 

We discuss the > case now in greater detail. Fig. 2 exhibits the allowed domain 
consistent with proton stability and relic density for the case mo = 700GeV as a function 
of At (where At is the t-quark Polonyi constant at the electro-weak scale) when tan /? 
and mi/2 are varied over the allowed range. Fig. (3) exhibits the allowed domain as a 
function of an (tan a// = cot (3) at At = when mo and mi/2 are varied over the allowed 
range of values. In each of the two cases one finds that the domain of the parameter space 
consistent with CDF, LEP data, proton stability and relic density is quite substantial 
even at the lower bound of 5 < 300 GeV"^ (M^g = 3Mg). 



New upper limits on the Higgs mass and on the chargino mass also emerge. One finds 
that nih < 105 GeV and m^^^^ < 100 GeV for B < 600 GeV-^ Thus the chargino should 
be seen at LEP2 while for much of the allowed parameter space, the light CP even Higgs 
should also be seen. The lightest ncutralino mass has an upper limit of 150 GeV and the 
maximum t-quark mass is ~ 165 GeV. These bounds are lower than the ones given in Ref. 
2 where no relic density constraint was imposed. The h, Wi and Zi mass bounds also 
decrease if one lowers the bound on B. 

IV. CONCLUSION : It is shown that in the physically interesting domain of the 
parameter space of the standard SU(5) supergravity, annihilation of the relic neutralinos 
is dominated by the light Higgs and the Z poles. Analysis is given which treats the 
thermal averaging over the poles rigorously. Previous approximate analyses are found to 
be inaccurate by several orders of magnitude near the Higgs pole and also significantly 
inacurate near the Z pole. It is found that for both fi > and fi < significant domains 
of the parameter space exists where all the desired constraints are satisfied. 
ACKNOWLEDGEMENTS : This research was supported in part by NSF Grant Nos. 
PHY-916593 and PHY-917809. One of us (P.N.) thanks TIFR where part of this work was 
done for the hospitality accorded him. We wish to thank Jorge Lopez and Kajia Yuan for 
pointing out an incorrect reheating factor (taken from Ellis et al., Ref. 14) in our earlier 
version of this paper. The correction leads to a much wider region where relic density 
constraints are satisfied. 

Figure Captions 

Fig. 1 : ^approx/^ as a function of mgiuino for top masses of 110 GeV (dashed curve), 
125 GeV (solid curve) and 140 GeV (dotted curve), showing massive breakdown of the 
approximation near the Higgs and Z poles. The poles occur close to where VLapprox/^ 
decreases sharply. Note that ^approx is least accurate in the region prior to the poles, 
which is also where Q.h'^ < 1. 

Fig. 2 : Allowed domains in the B — Af plane for top mass of 110 GeV (dashed curves), 
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125 GeV (solid curves) and 140 GeV (dotted curves) when mo = 700 GeV. The domain 
allowed by relic density constraints is the region between the upper and lower curves. The 
domain allowed by proton stability lies below the solid horizontal line when Mh^ < QMg- 
The gap in the central region for = 110 GeV is due to the requirement that rrih > 60 
GeV. 

Fig. 3 : Allowed domains in the B — an (tan an = ctn (3) plane for top quark masses 

of 110 GeV (dashed curves), 125 GeV (solid curves), 140 GeV (dotted curves) and 160 

GeV (dot-dash curves), when At = 0. The domain allowed by relic density constraints 

and proton stability is as in Fig. 2. 
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